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Aging in a two-dimensional Ising spin model with both ferromagnetic exchange and antifer- 
romagnetic dipolar interactions is established and investigated via Monte Carlo simulations. The 
behaviour of the autocorrelation function C{t, t^) is analyzed for different values of the temperature, 
the waiting time and the quotient 5 = Jo/ Jd, Jo and Jd being the strength of exchange and dipolar 
interactions respectively. Different behaviours are encountered for C{t,tw) at low temperatures as 5 
is varied. Our results show that, depending on the value of 5, the dynamics of this non-disordered 
model is consistent either with a slow domain dynamics characteristic of ferromagnets or with an 
activated scenario, like that proposed for spin glasses. 
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D ' Microscopic long ranged interactions are always of interest in different fields of physics because they can give rise 

^ to a variety of unusual macroscopic behaviours. Perhaps the better example in condensed matter are dipole-dipole 

^ I ^ ' interactions. In particular, the competition between long-range antiferromagnetic dipolar interactions and short-range 

' ferromagnetic exchange interactions can give rise to several interesting magnetic phenomena. Recent works in two 

■ dimensional uniaxial spin systems, where the spins are oriented perpendicular to the lattice and coupled with these 

^ ; ■ kind of iiitfiractions, have shown a very rich phenomcnological scenario concerning both its equilibrium statistical 

d . mechanicsu'B and non-equilibrium dynamical propcrticsa. Magnetization processes in these kind of systems are of 

^ ' interest due to aspects related to information storage in ultrathin ferromagnetic films. Moreover, there are several 

1 , contexts in which a short-ranged tendency to order is perturbated by a long-range frustrating interaction. Among 
' ^ others, model systems of this type has been pi'Ogftij>fd to study avoided phase transitions in supercooled liquids and 
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iJ = -Jo 5] a,a, + Jd^^ (1) 



^ \ charge density waves in doped antiferromagnetsHi 

^ ■ The above mentioned systems can be described by an Ising like Hamiltonian of the type 

'^^^J + 2^ -3 

■ where the spin variable Oi = ±1 is located at the site i of a square lattice, the sum X)<i j> runs over all pairs of 

, nearest neighbor sites and the sum j) runs over all distinct pair of sites of the lattice; r,y is the distance (in 

' crystal units) between sites i and j; Jo > and > are the ferromagnetic exchange and antiferromagnetic dipolar 

, coupling parameters respectively. For simplicity, we rewrite this Hamiltonian as follows: 

i' . with S — Jo/Jd- There are few numerical results concerning the equilibrium statistical mechanics, i.e. the finite 
' temperature phase diagram, of this model. In a recent work Maclsaac and coauthor^ have shown that the ground 
^ ' state of Hamiltonian (J^ is the antiferromagnetic state for S < 0.85. For S > 0.85 the antiferromagnetic state becomes 
O unstable with respect to the formation of striped domain structures, that is, to state configurations with spins aligned 
along a particular axis forming a ferromagnetic strip of constant width h, so that spins in adjacent strips are anti 
aligned, forming a super lattice in the direction perpendicular to the strips. They also showed that striped states of 
increasingly higher thickness h becomes more stable as S increases from 6 — 0.85. Moreover, they showed that the 
I striped states are also more stable than the ferromagnetic one for arbitrary large values of S, suggesting such a ph^g 
d ' to be the ground state of the model for 5 > 0.85. Monte Carlo calculations on finite lattices at low temperatureH'tl 
gave further support to this proposal, at least for intermediate values of S. Furthermore, such simulations have shown 
that striped phases of increasingly higher values of h may become thermodynamiGaily stable at finite temperatures 
for intermediate values of S. This results are in agreement with other analytic onesElEl. For low values of 5 the system 
presents an antiferromagnetic phase at low temperatures. At high temperatures, of course, the system becomes 
paramagnetic. 

The dynamics of the model is characterized by the formation and growth of magnetic domains, due to the compe- 
tition between the exchange and the dipolar interactions, which at low temperatures generate very large relaxation 
times. In a Monte Carlo study of the time evolution of the magnetization, Sampaio et aB have shown the existence of 
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two different type of relaxation, according to the value of 6. For 5 > Sc ^ 2.7 the magnetization relaxes exponentially, 
with a relaxation time that depends both on the temperature and 6. For 6 < 6c the magnetization presents a power 
law decay, with an exponent independent of 6. They also showed the existence of strong hysteresis effects in the 
presence of an external magnetic field for d > 6c- Since the relaxation time in this case increases with 6, the area of 
the hysteresis loops increases with 6. Hysteresis can also be detected in a simulation by looking at the evolution of 
the system starting from different initial configurations. In Fig.l we show snapshots of the evolution of a iV = 64 x 64 
system and 6 — 4:. For 6 < 6c hysteresis effects becomes very small and the magnetization presents a quasi linear 
dependence on the magnetic field, up to the saturation field (see Fig. 2). Hysteresis is a typical non-equilibrium phe- 
nomenon associated to domain dynamics with large relaxation times. In this scenario one also expects the presence 
of aging effects, that is, history-dependence in the time evolution of the response functions after the system has been 
quenched into some non-equilibrium state. |-| 

Aging phenomena has been vastly studied in disordered systems such as spin glasses (see Ref.tl and references 
therein), which are essentially out of equilibrium on experimental timpusealfis. However, they appear also in the phase 
ordering kinetics of ordered systems, such as the Ising ferromagnetBll3lll, associated with a slow domain dynamics. 
Aging can be observed in real systems through different experiments. A typical example is the zero-field-coolingE3 
experiment, in which the sample is cooled in zero field to a sub-critical temperature at time to- After a waiting time 
a small constant magnetic field is applied and subsequently the time evolution of the magnetization is recorded. 
It is then observed that the longer the waiting time tw the slower the relaxation. 

Although aging can be detected through several time-dependent quantities, a straightforward way to establish it in 
a numerical simulation is to calculate the spin autocorrelation function 



where < • • • > means an an average over different realizations of the thermal noise and is the waiting time, measured 
from some quenching time Iq = 0. 

In this work we present the results of Monte Carlo simulations in the two-dimensional Ising spin model defined by 
the Hamiltonian (|^) on a TV = 20 x 20 square lattice with free boundary conditions. We chose the heat-bath algorithm 
for the spin dynamics and time is measured in Monte Carlo steps per site. The quantity ^ is averaged over 100 
samples; for each run the system is initialized in a random initial configuration corresponding to a quenching from 
infinite temperature to the temperature T at which the simulation is done. We analyze the behaviour of C{t^t^) as 
a function of the observation time t, for different values of t^, , 5 and T. 

At enough high temperatures we find that the system does not present aging, that is, for any value of 6 there is a 
temperature above which C{t,tw) is independent of t^, as expected in a paramagnetic phase. At low temperatures 
we find different types of aging behaviours as 6 is varied. 

The typical behaviours of C{t,t^) are illustrated in Figures 3 and 4, for T — 0.5 and different values of 6 (waiting 
times tu, = 5"(n = 2, . . . , 6)). We also analyzed, for the same values of T and 6, the time evolution of the magnetization 
per site m(t) and the staggered magnetization per site ms{t) starting from different initial conditions , in order to 
characterize the different relaxation regimes. 

In Fig. 3a we see the typical behaviour of C{t,t^) in the antiferromagnetic state. The characteristic signature of 
C{t,t^) in this regime is the appearance of a plateau at some intermediate value of t independent of tyj [t ~ 10'^ for 
5 = —1), where C{t,tw) remains constant for a period of time that depends on T and 6; after such period C{t,tw) 
relaxes to zero. We also see a dependency on tw, that is, aging. The inset shows the evolution ofm{t) (filled circles) and 
ms{t) (open circles) starting from a random initial condition. We also analyzed the evolution of the same quantities 
starting from a fully magnetized initial state (to(0) = 1). We did not find any hysteresis effect in this region, the 
only difference being a much more slower convergence of ms(t) towards a constant value. These kinds of behaviours 
appear for negative values and also for small positive values of 6 at enough small temperatures, in agreement with 
previous calculations of the phase diagram of the modelcl. The existence of these plateaus suggests some type of 
quasiequilibrium state. This behaviour is rather unusual in this kind of system and it is probably related to the long 
range character of the interactions. 

In Fig. 3b we see the typical behaviour of C(<, tw) in the paramagnetic state, where it does not depend on tw In this 
case both m{t) and ms{t) converge fast to zero, independent of the initial configuration, as expected. This result also 
illustrates the characteristic shape of the low temperature phase diagram of the model (see Fig. 6 of RefE): following a 
line of constant temperature in the {6, T) diagram one can go from the antiferromagnetic region for low values of 6 to 
the striped one for high values of 6, passing through a paramagnetic region at some intermediate values of 5. For low 
enough temperatures it is expected a transition line between the antiferromagnetic and the striped regions, without 
an intermediate paramagnetic one. It is worth to note the logarithmic decay of C{t,tw)- This is a characteristic 
feature of this dynamics, that is, the relaxation of C{t,tw) is always slow, even in the paramagnetic region. 
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In Fig. 4 we show C{t,tw) for different values of S corresponding to the striped region, belowi (figures 4a and 4b) 
and above (figures 4c and 4d) Sc ~ 2.7. The behaviour of nis(t) in all these cases shows no traces of antiferromagnetic 
ordering. On the other hand, the behaviour of m{t) starting from a random (m(0) — 0) and a fully magnetized 
(to(0) = 1) states shows a clear differentiation between both regions: while for S > Sc stroiig hysteresis effects appear 
(see also Fig.l), for S < Sc such effects are negligible, in agreement with the results of Rcffl (see also Fig. 2). 

We found that C(t, i^) obeys a different type of dynamic scaling law for every one of the two striped regions, as 
can be seen from the data collapse of figure 5. For S < 6c C{t,tw) obeys the dynamic scaling 

C{t,t^) X cs {\og{t)/ \og{T{t^,))) , (4) 

(see figures 5a and 5b), while for S > Sc it obeys the following scaling 

C{t,U,)^cs{t/T{t^)), (5) 

(see figures 5c and 5d), where the time scale T(tu,)0 is shown in the insets. We observed, for values of S far enough 
of Sc, the scaling form T{tw) ^ i° . The log-log hnear fitting from Fig. 5 gives the following values: a — 0.83 ± 0.02 
for S = 2.0; a = 0.93 ± 0.05 for S = 3.5 and a = 1.00 ± 0.05 for S = 4.0. This values implies for S well below Sc a 
logarithmic scaling C[t,tw) oc c (log (t)/ log (t^,)) (Fig. 5a). The scahng form of C{t,tw) deviates from the previous 
one as S approaches Sc (Fig. 5b). For S > Sc we observe a scaling form C{t,tw) oc c{t/t'^) (figures 5c and 5d), with a 
exponent a that approaches the unity as S increases. 

It has been proposedLj'lliJ that aging phenomena are based on a slow domain growth at low temperatures, where after 
certain time t a characteristic domain size L{t) is reached. The time evolution of quantities like the autocorrelation 
function will then present a crossover from dynamical processes characterized by length scales smaller than the already 
achieved domain size to processes at larger scales dominated by domain growth through the movement of domain 
walls. In this scenario, scaling argumentsEJ leads to the following expected dependency of C{t,t^): 

C{t,t^)^c{L{t)/L{t^)), (6) 

at least for large values of t and tw 

Our results for S > Sc are consistent with a t/t^ scaling for large values of S, where the ferromagnetic short range 
interactions are dominant. This type of scaling, which is associated with an algebraic growth of the domain,. 
L{t) oc t^, appears in the slow domain growth dynamics of disordered systems with a ferromagnetic ground state 

Much more interesting is the logarithmic scaling form log (<)/ log (tj^) when S < Sc- This scaling is^associated with a 
logarithmic time dependence of the domain size L{t) oc (iogtw)'^ predicted by an activated scenariolj in spin glasses, 
in which disorder and frustration generate active droplets excitations with a broad energy distribution. 

In summary we have shown that the interplay between short- and long-range competitive interactions in an ordered 
system give rise to different kinds of aging at low temperatures. The different behaviours appear to be related to 
different domain dynamics as the relative strengths of the interactions are changed. In particular, for intermediate 
values of S, where the strengths are comparable, the competition between them generate a slow relaxation dynamic 
consistent with an activated scenario, like that proposed for spin glasses. This result suggests a possible relationship 
between the microscopic dynamic properties of this model and that of disordered frustrated systems. In this sense, 
it would be of interest to investigate, for instance, the possible existence of broad energy distributions of low-lying 
excitations in the present model. Works along these directions are in progress and would be published elsewhere. 
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FIG. 1. Snapshots of the time evolution of a 64 x 64 system with 5 = A, T = 0.5 for different initial conditions. Figures (a), 
(b) and (c) correspond to an initial configuration with mean magnetization per site m(0) = 0, for t = 1, 10 and 1000 Monte 
Carlo steps per site respectively. Figures (d), (e) and (f) correspond to an initial configuration with m(0) = 0.6 and the same 
values of t. 



FIG. 2. Snapshots of the time evolution of a 64 x 64 system with S — 2, T = 0.5 for different initial conditions. Figures (a), 
(b) and (c) correspond to an initial configuration with mean magnetization per site m(0) = 0, for t — 1, 10 and 1000 Monte 
Carlo steps per site respectively. Figures (d), (e) and (f) correspond to an initial configuration with m(0) = 0.6 and the same 
values of t. 



FIG. 3. Autocorrelation function C{t,tm) vs the observation time t at T = 0.5, for = 5^ (circles), = 5^ (squares), 
tw = 5* (triangles) and t„, — 5^ (hexagons), (a) 5 — —1, inset: magnetization m{t) (filled circles) and staggered magnetization 
nisit) vs time, (b) S = 0.5 



FIG. 4. Autocorrelation function C{t, t^) vs the observation time t for — 5^ (circles), = 5* (squares), = 5^ (triangles) 
and tw = 5® (hexagons), for different values of S corresponding to the striped phase. Insets: magnetization m{t) vs time starting 
from random (m(0) = 0, open circles) and a fully magnetized (m(0) = 1, filled circles) states. 



FIG. 5. Data collapse of the curves shown in Fig. 4 (the values of and the corresponding symbols are the same in both 
figures) . 
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Figure 5 



